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Abstract 

Given a regular polygonal arrangement of identical objects, turning 
around a central object (masses, vortices or dNLS oscillators), this pa- 
per studies the global bifurcation of relative equilibria in function of a 
natural parameter (central mass, central circulation or amplitude of the 
oscillation). The symmetries of the problem are used in order to find 
the irreducible representations, the linearization and, with the help of a 
degree theory, the symmetries of the bifurcated solutions. 
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1 Introduction 

Consider a polygonal arrangement of n identical objects turning in a plane, at 
constant angular speed, around a central object. These objects may be masses, 
following Newton's law of attraction, or point vortices, with Kirchoff's law, or 
nonlinear oscillators coupled to nearest neighbors in a finite circular lattice and 
a common phase. 

A relative equilibrium for these problems is a stationary solution in rotating 
coordinates, [9]. For each angular speed there is a regular polygonal relative 
equilibrium and an associated central quantity (mass, or circulation or ampli- 
tude of the oscillation) which is taken as a parameter. 

The purpose of this paper is to prove that for certain explicit values of this 
parameter there is a global bifurcation of relative equilibria with a specific sym- 
metry. The tools for this study is representation theory and a simple version 
of the equivariant topological degree, studied in j7|. As a matter of fact, the 
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reduction to irreducible representations gives a very clear picture of the sym- 
metries involved and will be used in forthcoming papers on the bifurcation of 
periodic solutions for these problems. 

There is a vast literature on the (n + l)-body problem, beginning with 
Maxwell's model for the Saturn rings. The point vortices problem has also 
attracted a good deal of research. This is usually done with a combination 
of numerical and explicit computations, where the symmetry is regarded as a 
nuisance. In the present paper we hope to show that these symmetries, when 
considered as a whole, facilitate instead the study. 

In the rest of this introduction we shall set more precisely the problems. 
Then, in Section two, we shall see how the symmetry of the problem forces 
the Hessian of the system of equations to have a special structure and we shall 
introduce a transformation which will bring the Hessian in a block-diagonal 
form, according to the different isotropy types. Afterwards, in Section three, we 
shall state our bifurcation results, both local and global, giving solutions with 
specific symmetries. In Section four, we shall complete the spectral analysis for 
the n-body and n-vortex problems. The final section is on the discrete NLS, 
which belongs to a somewhat different field of application but where a very 
similar analysis may be performed. 

Among the papers listed in the bibliography, in particular [TU], [13], [T], 
[TT] and [12] and their respective references, the paper closest to our results is 
[TU] for the n-body and the n-vortex problems. These authors find the same 
critical values of the parameters and use a normal form analysis and numerical 
computations in order to prove local bifurcation results. 



1.1 {n + l)-vortex problem 

One of our purposes is to study relative equilibria of n -I- 1 vortices with circula- 
tions fj,o = fi and fj,j = 1 for j G {1, n}. Let qj{t) € be the position of the 
vortices and Xj{t) — e~'^'^*qj{t) be their position in rotating coordinates, with 
angular speed uj. Then, the dimensionless equations for relative equilibria with 
frequency w are 

n 

EXj Xi 

where J is the canonical symplectic matrix. 



1.2 {n+ l)-body problem 

Another of our purposes is to study relative equilibria of n + 1 bodies in the 
plane, where the bodies have masses fiQ = ii and fXj — 1 for j e {1, ..,n}. Let 
Xj(t) be the position of the bodies in the plane and in rotating coordinates, 
as above with angular speed (w)^/^. It is well known that, after the change 
of variables, the dimensionless equations for relative equilibria with frequency 
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(a;)^/^ are 

n 

1.3 General problem 

Now, we will set a formulation generalizing both previous problems. Let x be 

the vector (xojXi, ...,x„)-^, where T denotes the transposed, and Ai the matrix 
diag{iJ,, 1, 1). Our aim is to look for critical points of the potential 

V«(x) = uj{x'^Mx)/2 + ^iiiiJ,j(j)o,{\\xj - Xi\\), (1) 

where 4'a{x) satisfies ^^(x) = —l/x°' for a € [1, oo). 
Since the potential V has gradient 

n 

VxjV^x) = ujiJjXj - HiHj 



i=0 II^J 

then the critical points of V are the relative equilibria of the vortex problem 

for a = 1, and of the body problem for a = 2. Also, the case a > 1 can 
be regarded as a problem of relative equilibria for bodies, where the general 
attraction potential is <^oi- 

Hereafter we represent points in and C indistinctly. Let C, =2TT/n and let 
us set the positions of the bodies aj as ao = and aj = e*^'» for j € {1, n}. We 
see that aj form a relative equilibrium with a central massive body at the origin 
surrounded by bodies of equal masses in a regular polygon. This polygonal 
relative equilibrium was studied by Maxwell as a simplified model of Saturn 
and its rings. 

Proposition 1 a = (ao,...,an) is a critical point of the potential V{x), when 
w = /U + si with 



_ ^ 1 - e^J"- _ 1 V- 1 

||l-e'jqr+' ~ 2^ t:^ sin("-i)(i, 



,=iir---|i - ,=1— (iC/2) 

Proof. For j = 0, we have Vxo^(^) = A* Y^]=o ^^^^ — 0- 3 0? have 

n 

Vxj V{a) = ujaj - ^ - naj = aj (w - (m + si)) , 

Therefore, a is a relative equilibrium for frequencies to — /i + Si. ■ 

Notice that any homograph of a relative equilibrium is also a relative equi- 
librium. This is why we have decided to fix the norm of the relative equilibrium 
and leave the parameter jj free. Our objective is to find global bifurcation of 
relative equilibria from aj using the parameter fj,. Now let us see the symmetries 
of the problem. 
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Definition 2 Let Sn be the group of permutations o/{l,...,n} and let Dn be 
the subgroup generated by the permutations (^{j) = j + 1 and K,{j) = n — j. We 
define the action of Sn in M2("+^) as 

p{-y){xQ,Xi.....Xn) = (xo,X^(i),...,X^(„)). 

In addition, we define the action of 0{2) = U kS^ as 

p{0) = e~^'^ and p{k) = U, 

where TZ is the matrix diag{R, R) with R = diag{l,—l). 

Because n of the bodies have equal masses, the potential V is 5„-invariant. 
Moreover, the potential V is 0(2)-invariant since the equations are invariant by 
rotating or reflecting the positions of all the bodies. Consequently, the gradient 
W is F-equivariant with F = 5„ x 0(2). This means just that 

WV{p{j)x) = p{j)VV{x) 

for all 7 e F. 

Let Dn be the group generated by the elements {(, Q and (k, k) of Sn x 0{2), 
where C = ^n/n e 5^ The action of (C, C) and (k, k) in ^^(n+i) jg 

(C, C)^; = p{0^~^^x and (k, k)x = p{k)TZx. 

As the action of (C, C.) and (k, k) leaves the equilibrium a fixed, then its isotropy 
group, i.e. the subgroup of F which fixes the orbit a, is 

r-a = Dn. 

2 Irreducible representations 

In order to prove the bifurcation theorem, we need to find the spaces of irre- 
ducible representations of 

Let us define Aij to be the 2x2 submatrices of D'^V{a) such that 

D'V{a) = A={A,j)2j^,. 

Due to the fact that D'^V{a) is i?„-equivariant, one has the following result: 

Proposition 3 The blocks Aij satisfy the relations 

Aij = e""''^A^(i)f(j)e'''^ and Aij = RA^(^)^(^j)R. (2) 

Proof. Since the matrix D'^V{a) is £)„-equivariant, then the matrix A and 

p{C,)e.~'^^ commute. Therefore, 

A = p{Qe-^iAe^^p{Q-\ 
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Hereafter, we denote by [u]i the coordinate Ui G of the vector u = {uq, 

Therefore, 

[p{Oe--^<Ae-^ip{0-'u]i = e-'<[Ae'^<p{C)-'u]^^,-, 

From this equaUty, we get that 

J2AijUj = [Au]i = ^e^^^A(;(,)^(j)e''fuj-. 
j 

Then we conclude that Aij = e~'^^A^(j')|j(-j')e'^'». Using a similar argument and 
the fact that A and p{k)TZ commute, we obtain Aij = RA^f^i^^f^j^R. ■ 

Now, we may find the irreducible representations of the action of the group 

Dn. 

Since the irreducible representations are different for n = 2 and n > 3, we 
shall concentrate of the case n > 3 in the remaining part of the paper, except 
for comments on the case n = 2. 

Definition 4 Let us define the vectors vi and Vn-i as 

For k e {2, n — 2, n}, we define the isomorphisms : ^ Vk as 
Tk{z) = (0,n-i/2e(*'=^+^)^z,...,n-i/2e"(*'=^+^)C^) u,ith 
Vk = {(0, e(*'=^+^)^2, e"(''=^+-^)f2) : z G C^}, 
and for k € {1, n — 1\ , we define the isomorphism : C'^ — > Vk as 
Tk{a,w) = (wfca,n-i/2g(<fe/+./)Cy;^_ _^„-i/2g«(»fe/+j)Cy,) ^^^^^ 
Vk = {Ka,e(''=^+^)^u;,...,e"(^'=^+'')«u;) :w€C^,a€ R}. 
Next let us find the action of the group Dn on the subspaces Vk- 
Proposition 5 The actions of {(, () and (k, k) on Vk are 

{C,OTkiz)^Tk{e'''^z) and 

{K,K)Tk{z)=Tn-k{Rz), 

where R is the matrix diag{l, 1, —1) for the special cases A; e {1, n — 1} . 
Proof. For k G {2, .., n — 2, n}, we have 

= [Tk{e^^'''+'^'^z)]j. 
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Therefore p{C)Tk{z) = Tk{e^'^^+-^'^'^ z). Since the element ( G 0(2) acts as 
e-^^Tk{z) = Tk{e-^-'z), we conclude that (C, C) acts as (C, OTk{z) = Tk{e''"^z). 
For k e {l,n — 1} we have, as before, that 

Moreover, from the equality e~"^'»Vfc = e^'"'Vk we get that ( G 0(2) acts as 

e-^-^Tk{a,w) =Tk{e"'ivka,e-'^'^w). 

Hence, the action in this case is also {(,C,)Tk{z) = Tk{<'^'^^z). 

It remains to find the action of (k, k). For k e {2, n — 2, n}, we have 

[{K,K)Tk{z)]j = [7^Tfe(^)],(,) =n-i/2g,(-ifc/+j)c^^ 

therefore the action is (k. K)Tk(z) = T„_fc(i?z). For k G {1, n — 1}, by a similar 
argument and the fact that Rvk = ^^(fe), we prove that the action is as before 
but with i? = 1, —1). ■ 

Consequently, we have that the spaces Vk ® Vn-k are subrepresentations of 
the action of £>„. Moreover, the action of (C,C) and {k,k) on the subspace 
Vk e Vn-k is 

{C,C)izk,Zn-k) = {e'^^Zk^C'^^Zn-k) and 

{K,K){Zk.,Zn-k) = {RZn-k,RZk). 

Let Z„ be the group generated by (C,C)- Since the subspaces Vk are irre- 
ducible representations of Z„, by Schur's lemma (that is a linear map which 
commutes with action, must send equivalent representations into themselves), 
we obtain D^V{a)Tk{z) = Tk{Bkz). Furthermore, as D'^V{a) commutes with 
the action of {k,k), then the blocks Bk must satisfy BkR = RBn-k- Conse- 
quently, there must be a map that puts the matrix D'^V{a) in diagonal form with 
the blocks Bk- Clearly, the isomorphisms Tk, with range Vk, are the components 
of this orthogonal transformation. 

Proposition 6 Define the map Pz = X)fe=i Tk{zk) for z = (zi, z„), then the 
linear map P is orthogonal P* = P~^ . 

Proof. Since the matrix e^^^ is an isometry in and 

Y^e^^i^-^)i=n6ki, 

3=0 

for k, I G {2, n — 2, n}, then 

n 

{Tk{zk),Ti{zi)) =n-'J2^''^''~''"^ {e''^zk,e^'<zi) = 5ki {zk,zi). 
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In fact, since vi and are orthonormal vectors, one proves that {Tk{zk), Ti{zi)) 
6ki {zk,zi) for all k and I. Therefore, the map P satisfies 

Thus, P is an isometry on C2("+i) and P* = p-\ m 

By Schur's lemma, the matrix D^V{d) is diagonal in the new coordinates, 
that is 

p-^D^V{a)P = diag{Bi, B„). 



Om' next objective consists in finding the blocks Bk in terms of the matrices 
Aij . Remember that the matrices Aij are the 2x2 submatrices of the Hessian 
D^V{a). 

Proposition 7 For fce{2,...,n — 2,n} the blocks Bk are 

n 

Proof. For ; ^ we have [ATk{z)]i = nr'^^^ Y."=i ^ij^^''''''^''^^ Now, from 
the relation ([2]), we prove that Aij — e''^^A„(j_;')e^'"'^, with I — j modulo n. 
Hence 

n 

Consequently, we rewrite the sum as 



[ATk{z)]i = 



- „-i/2J(ife+./)C 



gj(ifc/+,/)C^ 



= [Tk{Bkz)]i 



But since the isomorphisms T^, are defined on the invariant subspaces 14, then 
[ATk{z)]o = [Tk{Bkz)]o and we conclude that ATk{z) = Tk{Bkz). Actually, one 
may prove directly that [ATfe(z)]o = [Tk{Bkz)]o, for instance see [1]. ■ 

Proposition 8 For k E {l,n — 1} the blocks Bk are 



Bk 



^ l^/^Ar^^Vk Y.%iAr.,^'^^^'^'^<^ )' 

Proof. For / 7^ 0, we have 

n 

[ATk{a,w)]i = (A,o«fc)a + n-i/2^^i,e-''(^'=^+^)««;, 

where a £ C and w G C'^. Now, from the symmetries ([2]) we prove that Aiq 
e'-^^Anoe'^'^''. Hence, Ai^Vk — e''^^yl„oe"''^^Wfe. Moreover, since e~'''''Wfc 
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then AiQVk = e''^'^*"'"'^-'^A„o'yfc- Using the previous computation, we find 



that 

[ATk{a,w)]i = n 



^ „-l/2gi(«fc/+,7)C 



(3) 



For / = 0, we have [ATk{z)]o - (AooWfc) a+n-i/2^^y,^ ^ith = Ej=i ^0^6^^*'=^+''^^. 



From the relations ©, we have that Aqq = cl. Now, since Vkv'J\ = 1, then 



[ATk{a,w)]o ^ Vk (efAooei)a + n ^^^{v'^Dk)w 



(4) 



Consequently, from the equalities ([3]) and (jl]), we obtain Z)^T^(a)Tfc(a, w) 
Tk{Bk{a,w)), with 



n^'^A^^Vk EU An,e^^''''+''>^ 



Moreover, since the map P is orthonormal and the matrix A is selfadjoint, 

T 

then Bk must be selfadjoint and {v^Dk) — nAnoVk- Actually, one may prove 

T 

directly that nAnoVk — (vjDk) ~ ^k'^-'k, for instance see [4]. ■ 

Remark 9 In the computation of the blocks Bk, we have used only the sym- 
metries of D^V{a). This will enable us to apply these results to a wide class 
of problems, as the dNLS equations at the final section. Also, notice that the 
change of variables was done in complex coordinates, and these will allow us to 
prove bifurcation of periodic solutions in a series of forthcoming papers analo- 
gous to J^: as a matter of fact, the natural approach to the study of periodic 
solutions is, in this context, the use of Fourier series. Thus, the change of 
variables, which we have introduced, will be helpful. 

However, in order to find bifurcation of relative equilibria, we need the 
change of variables for real coordinates. 

Proposition 10 // the matrix D^V{a) has domain ]R.^("+-'^', then the matrix 
p-^D^V{a)P has domain W = p-i]R2("+i) and 

p-^D^V{a)P = diag (Si, B2, B„/2, B„) with 
W = X X ... X X m2. 

Moreover, the action on the block Bk is 

(C, C)^*; = e^'^'Zk and (k, K)zk = Rzk, 

where R is the matrix diag{l,l,~l) when k = 1 and diag(l, —1) for the re- 
maining cases. 
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Proof. First, we need to identify the subspace W = {z : Pz G M2(n+i)| jf 
is real, then 

n n 

Y,Tk{zk) = Pz = T^ = Y,Tn-k{Zk). 
fe=l fc=l 

Thus, the subspace W is the set of points (zi, Zn) such that z„_fc = Zk- 

Remember that (k, k) acts on the coordinate Zk as (k, K)zfc = Rzn~k- Hence, 
for k € {n/2,n} we have Zk = Zn-k = Zk & K^, then {k,k) acts as {K,K)zk = 
Rzk- Consequently, the blocks -B„/2 and Bn are defined in a real space with 
real action. 

Now for k ^ {n/2, n}, we can take the isomorphism T{zk) = {zk, Zn-k) with 
Zn-k = Zk- In this way [k, k) acts as 

{K,K)T{zk) = {Rzk,Rzk) = T{Rzk). 

Moreover, since ((, () acts as (C, C)£fe = e^'^^Zk, then {(, C)T{zk) = T{e^^^Zk). Fi- 
nally, we use the equality Bn-k = Bk to prove that {Bk, Bn-k)T{zk) = T{BkZk)- 
■ 

Remark 11 If n = 2 we have to define T2 as before, from into V2. 
However, for k = 1, define the isomorphism Ti : ^ Vi as 

Ti{v,w) = (w,2-^/^u;,2-^/^u;) con 
Vi = {{v, 10,10) : v,w e C^}. 

Then, tha action of D2 on Zk &Vk is 

(C, 0^2 = Z2 y {k, k)Z2 = RZ2, 

(C>C)-2i = -zi y {k,k)zi = diag{R,R)zi. 

Hence, the spaces Vk are irreducible for the action o/(C, C), but Vi contains 
two representations, one where (k, k) acts as the identity and the other where 
this element acts minus the identity. 

One m,ay prove that Aij are diagonal matrices ans satisfy An = A22, A21 = 
A12 and Aqi = Aiq = ^402 = ^20- Thus, 

^00 = (si + lJ.)lJ.I - 2A20, A20 = ~iJ,diag{a, -1), 

^21 = --^^diag{a, -1) y A22 = (si + /J,)/ - (A20 + ^21)- 

In particular, contrary to the case n > 3, where Aqq is a multiple of the 
identity, this matrix is only diagonal. 

The transformation P is orthogonal and one sends the hessian into diag{Bi, B2), 
with B2 = {a + + si)diag{l, 0), but Bi is now 

/ IJ,{si+ iJ, + 2a) -V2aiJ, \ 

A< (si + M - 2) V21J, 

-\f2aii si + (a + l)/i 

V \f2\x si / 
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In order to study the bifurcation of relative equilibria we need to restrict 
the block Bi to the fixed point suhspace of {k,k), that is to the first and third 
coordonates. There, 




( 



fi{si + fi + 2a) 



si + 



(a + l)M ) 



The determinant of this matrix is /i (/i + si) (2a + + si + a/i). At^ = 
0, —Si, there will be a bifurcation, as in the case n > 3, and at 



there will be a bifurcation with symmetry Di. For the vortices, one has a — I, 
Si ^ 1/2 and /ii = —5/4, while, for the masses, one has a = 2, si — l/A and 
^ll = -17/12. 

3 Bifurcation theorem 

We shall now give sufficient conditions for the bifurcation of relative equilibria 
from a. To carry on this proof, we need to apply the change of variables P 
directly for the potential V. 

In this manner, we define the potential Vp : — )■ R as Vp{x) = V{Px). 
Then the potential Vp is F-invariant and the gradient VVp is F-equivariant with 
the action pp{'i) = P^^p{'j)P. Note that xq = P^^a is the relative equilibrium 
in the new coordinates. 

To find the symmetries, for each h dividing n, we define the group Dh as 
the one generated by the elements {n/h){(, () and (k, k). These groups Dh are 
subgroups of the isotropy group £)„. Our approach consists in applying Brouwer 
degree to the maps \i'Vp{x) restricted to the spaces of fixed points of Dh- As 
seen in [7], this is equivalent to the Dn-equi variant degree. 

So we set the function fh{x) as 



Then, the zeros of fh{x) are the relative equilibria with symmetry Dh- Now, the 
polygonal relative equilibrium is xq — P^^a, so xq is a zero of fh{x)- Since we 
wish to prove existence of bifurcation from xq, we need the sign of det /^(xq). 

Proposition 12 Define Ck as 



fii ^ - {2a + si) /{a + 1) 



fh{x) = VVp{x)\^n, --W^-^W 



(Tfe — sgn (e^-BfeCi) for k G {n,n/2} and 
cTfc = sgn (det Bk) for A: G [1, n/2) n N. 



(5) 



Then 



sgn {detfh^xo)) 



= cr, 



n 



n 



a 



je[l,ri/2]nNh 
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Proof. Since (C, C) acts on the coordinate as e^^^Zk, the action of {n/h){(, C) 
on Zk is trivial when This happens precisely for the coordinates 

Zk with k € hN. Now, the action of (k, k) on the coordinate is trivial whenever 

Zk — Rzk- Therefore, z S W^*^ only when z — {zh, Z2h, Zn) with Zk — Rzk- 
Then, the matrix D'^Vp{xa), on the space W^'', is 

diag{Dh,D2h,-;Dn), 

where the blocks Dh are as follows: 

• For k e {n,n/2}, we have that Dk = ejBkei because z/j e R x {0}. 

• For k ^ {1, n/2, n}, since z^ S R x iR, we have that Dk ~ T*BkT, where 
T — diag{l, i) is the natural isomorphism between R^ and R x iM.. 

• For fc = 1, since z/, G x iR, we have Dk = T*BkT, where T = 
diag{l, 1, i) is the natural isomorphism between R'^ and R'^ x iR. 

Finally, from the definition of <7k, we get that sgn (detDk) = crk and there- 
fore 

sgn det fp{xo) — sgn detZ?„ dctl?j = 

\ je[l,n/2]nNh / 

■ 

Hence, we have given the sign of det flixo) in terms of the blocks Bk- 
3.1 Local bifurcation 

In order to apply Brouwer degree and prove bifurcation, let us define /, from 
B2e X to R X VF, as 

f{x,fi) = {\\x - xo\\ - ejh{x,fi)) with 

B2e X B2p = {(x, ^I) e W^'^ xR:\\x-xo\\<2e,\n-fio\< 2p}. 

Theorem 13 The Brouwer degree of f{x,ii) is well defined and 

deg(/; B2e x B2p) = ?7h(Aio) = nhipo - p) - rihipo + p). 

Hence, when r]h{po) 7^ 0, there is a local bifurcation from (xq, /io) with symmetry 

Proof. As in [6 , the proof consists in a linear deformation of the function 
||a; — xoll — e to p — 11^ — poll and of the function fh{x) to — xq). Then, we 
may use the excision property to prove that 

deg(/(p);B2pXS2e) = deg{f'i^{po-p){x-xo); B2e)-deg{f'f^{po+p){x-xo); B2e). 
Therefore, the first part of the proof follows from the fact that 
deg(/;l(p)(x - xo); B2e) = n/i(p). 
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Now, supposing r/hino) 7^ 0, for small e there is a (x^,^^) with C W^'^ 
such that fh{xe,fJ:s) = and d{xe,xo) = s. Moreover, when we let s tend to 
zero, by the compactness we have a series £fe ^ such that ^ AH- By the 
continuity we conclude that /ii = /iQ. ■ 

When only one of the blocks has a determinant which changes sign, we 
have the following result. 

Theorem 14 For k G {1, [n/2], n}, let h be the maximum common divisor 
of k and n. Supposing ak{fi) changes sign at no and (Jjino) /or the others 
j € [1, n/2] n (N/i), then there is a bifurcation with maximal symmetry D^. This 
means that the local bifurcation is in U^^^^ . 

Proof. To assure all the symmetries of the bifurcation, we apply the previous 
theorem with h the maximum common divisor of k and n. By hypothesis the 
product 

je[i,n/2\rmh {jjtk) 

is not zero, then r?/i(/Uo) = ±2. Henceforth, there is a bifurcation in the fixed 
point space of Dh- 

It only remains to prove that Dh is the maximum group of symmetries. Let 
Dp be a group such that Dh C Dp, this means that h divides p. Since p does 
not divide k, then 

sgn det fp{xo) = Oj 7^ 0. 

iG[l,n/2]nNp 

Consequently, the linear map /^(xo) is invertible, and by the implicit function 
theorem, we deduce the non existence of solutions in W^^ near (xq, /xq). ■ 

3.2 Global bifurcation 

Now, we wish to prove a global bifurcation result, which is just an adaptation 
of the Rabinowitz' alternative. Notice that this approach may not give all 
the best information available for a global result, as an application of the F- 
cquivariant degree. But this equivariant degree (for this larger group) presents 
strong technical difficulties. 

Let us define T as the set {{Txq,\x) : e M} with xq = P~^a. Let S be 
the set of zeros of fhix,^), wc define G = S\T as the nontrivial solution set. 
G\G C T is the bifurcation set and an element of (xq, /xq) S G\G is said to be 
a bifurcation point. Let C C (5 be the connected component of the bifurcation 
point (xo, ^o). Then C CiT consists of the bifurcation points of the branch C. 

We define the collision set as 

4- = {x e . ^. ^ ^. ^ j-)!^ 

Also let Ap = {||/x|| < p} and ftp be 

= {x e M2(n+i) . ii^ii < p^p-i < d{x,^)}. 
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Since, for p big enough, the set Q.p is a big ball without a small neighborhood 
of hyperplanes of codimension 2, then the set J7p is connected. 

We say that the component C is admissible whenever it is contained in some 
set rip X Ap. Otherwise we say that C is inadmissible and this corresponds to 
the cases where (a): the parameter fi, on the component, goes to infinity, or (b): 
the norm of x on the component goes to infinity or (c) : the component ends at 
a collision point. 

Theorem 15 // the component C is admissible and the set C HT is isolated, 
then C returns to other bifurcation points {xr,Pr)} o,nd 

r]hixo,p.Q) + ■■■+r]h{Xr,Pr) ^0- (6) 

Proof. Since C is admissible, we may construct a set C rip x Ap such 
that C C with dCl O C = (/) and such that fh{x,^) is zero on 9f2 only 
when X Cz T. Since fh{x,fj,) is not zero on dfl unless x £ T, then the degree 
deg{d{x, T) — e, fh\ ^) is well defined. Moreover, as f2 is bounded, we can take 
£ big enough in such a way that this degree is zero. 

By hypothesis CRT consists of isolated points. Consequently, taking e small 
enough, the points which satisfy d(x,T) — e and fh{x,p) = in are in the 
finite and disjoint union of i?2e(a^o) x ^2p(Mo) for {xa,po) E C C\T. Hence, by 
the excision property of the degree we have 

= deg(d(x,r) - e,/ft;rj) = ^ deg(d(x, xq) - £, A; x Bap). 

(a:o,Mo)eCTlT 

We conclude the result from the computation of the local degree. ■ 
3.3 Symmetries 

The relative equilibria with symmetry D}^ are composed of bodies arranged as 
n/h regular polygons of h sides with some polygons related by reflection. To 
give a sharper description, let us call an h-gon as the set of positions 

{re''^e'=(2"/'') : yfc = l,...,/i}, (7) 

and a 2/i-gon as 

Proposition 16 In a relative equilibrium with symmetries Dh the central body 
stays on the real axis if h = 1 and remains at the origin if h > 1 . The other 
bodies satisfy the following arrangements: 

(a) If n/h is odd, the relative equilibrium has an h-gon ^ of bodies with p = 0. 

The remaining bodies form 2h-gons with Lp £ (0,7r//i). 

(b) If n/h is even, the relative equilibrium has two h-gons ^ of bodies, one 

with p = and another with ip — n/h. The remaining bodies form 2h-gons 
(0), with ip e (0,7r//i). 
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Proof. Since the central body satisfies the symmetry xq = {k, k)xo = xq, then 
xo € M. Moreover, as xo = (n/ft,)(C, C)xo = e~*^^/'*a::o, then xq = whenever 
h>l. 

Now, for the remaining bodies, j € 1, ...,n, we use the notation Xj — Xj^kn. 
Hence these bodies satisfy the relations 

(a) Xj ~ [n/h){QX)xj — e~^^'^^^^'^Xjj^n/h and (b) Xj = {K,K)xj = x^j. (9) 

By (a), the positions of the bodies are determined only by the bodies with 
j € Nn{—n/2h,n/2h], and by (b), these are determined by the bodies j G 
Nn[0, (n/2/i)]. 

From (a) we have an ft,-gon ([7]), for each j £ {0, (n/2h)}. Actually, from (b), 
we deduce that (p = for j = and (p = ir/h for j = n/2h. Now, for each 
body, j € Nn(0, (n/2/i)), we have a 2/i-gon Furthermore, since a 2/i-gon 
has collisions for ip G {0, ir/h}, then we can chose (p G (0, ir/h). ■ 

In order to give examples of the previous descriptions, we shall analyze the 
cases Di, D2 for n even, and for n = 6. 



I ---- 

! ■ > ■ — r 



(a) n = 3. (b) n = 5. 

Figure 1: Symmetries of the group Di. 

The group Di is a subgroup of and it is generated by (k, k). The relative 
equilibria with symmetry Di have the central body on the real axis. The other 
bodies satisfy the following: 

(a) If n is odd. One body is on the real axis. The remaining bodies form 

symmetric couples with respect to the real axis, (see the examples n = 3 
and n = 5). 

(b) If n is even. Two bodies are on the real axis without any relation. The 

remaining bodies form symmetric couples with respect to the real axis, 
(see the examples n ^ A and n = 6). 

The group D2 is generated by (tt, tt) and (k, k), and it is a subgroup of Dn 
whenever n is even. A relative equilibrium with symmetry D2 has the central 
body standing still at the origin. The other bodies satisfy the following: 
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(a) Symmetries of Di. (b) Symmetries of D2. 

Figure 2: n — 4. 



(a) If n/2 is odd. One pair of bodies is on the real axis symmetric with respect 

to the imaginary axis . The remaining bodies form squares symmetric 
with respect to both axes, (see the example n — 6). 

(b) If n/2 is even. One pair of bodies is on the real axis symmetric with 

respect to the imaginary axis. Another pair of bodies is on the imaginary 
axis symmetric with respect to the real axis. The remaining bodies form 
squares which are symmetric with respect to both axes, (see the example 
n = 4). 




(a) Symmetries of Di. (b) Symmetries of D2. (c) Symmetries of D3. 

Figure 3: n = 6. 

Finally, for the subgroup D3 of Dq we have n/h = 2. Therefore, the central 
body stands still at the origin and the remaining bodies form two triangles 
without relation, one with = and the other with (p = 7r/3. 

Note that these figures are only an illustration of the possible configurations 
which may happen. They have to be taken in this perspective, as it the case 
of the figures in other papers, such as [10]. A precise numerical analysis of the 
positions of the bodies far from the relative equilibria is outside our present 
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concern. Furthermore, our last proposition is mathematically valid for any /i. 
However, for the gravitational problem, the masses need to be positive or, at 
least if one is considering an attraction given by charges, that w should be 
positive, since we took its square root. Finally, the spectral analysis and some 
of the following remarks will give a complement of information and a better 
justification of our figures. 



4 Spectral analysis 

It is time to calculate explicitly the bifurcation points for the general potential 
(dl . We begin by computing the matrices Aij . 

Proposition 17 Define a+ = (a + l)/2 and a_ = (a — l)/2, then, for n> 3, 
we have 

Aqq = /i (si + a* + a-n) I, 
■^no — —fJ-{a_I + a_|_i?) and 

n-l 

Ann = {Sl + M)^ - X! ^"-J - 

In addition, we have for — 1} that 

(2sin(jC/2))"+'^ ' 
Proof. Notice that Wxi4'{\\xi — Xj\\) = —'^Xj4'{\\^i ^ Xj\\) for i ^ j, thus, 

Aij = fJ.ttJ,jDx^'^x,H\\0'i - flill) = -t^iy']Dl^cj){\\ai - aj\\). 
And, for the matrix An, we have 

Au = (si + + '^^ii^.jDl^cf){\\ai ~a.j\\) = (si + -'^Aij. 

Let us set aj = {xj,yj) and = \\{xi,yi) — {xj,yj)\\, then the function 
<^Q (dij ) has its matrix of second derivatives 

^ ^.[d,,) - y _ ^^^^^^ _ _ ^^^2 j ^,1- 

Since the distance from oq ~ (0, 0) to a„ = (1, 0) is dno — 1, then 



.0 -1 . 

Moreover, as e^^"'^ = for C 7^ tt, or equivalently n > 3, then 

n n 
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Therefore 

n 

Aoa = (si + ^J■)^J■I - ^ Aoj = ^ (si + /i + a_n) /. 

It remains only to find the matrix Anj for j e {1, n — 1}. As a„ = (Ij 0) 
and flj = (cosj'C, sin j'^'), then the distance dnj satisfies 

dl^ = (1 - cosjC)' + sin^ jC = 4sin2(jC/2). 

Using the previous results, we have 

a + l f (1-cosjC)^ -(1 - cosjOsinjC \ 
" d'^+' [-{^- cos jC) sin jC (sin jC)' J + <+i " 

Now, since sin^ jC = (1 — cosjC)(l + cosj^) and = 2(1 — cosjC), then 

_ 1_ / a + l f 1-cosjC -sinjC 



V 2 V -sinjC 1 + cosjC 



Finally, using dnj = 2sin(jC/2) we conclude the result. ■ 

It can be seen that the bifurcation points are just the points where det Bk 
changes sign. Now we can find explicitly the blocks Bk for the general potential 

O- 

Proposition 18 Define Sk, oik, o-nd 7^ as 

^ 1 sin^(fcjC/2) 
2"^4^sin"+i(jC/2)' 

ak = —{sk+i + Sk-i), Pk = a+isk - si) and jk = -^(sfc+i - Sk-i). 
For k e {2, n — 2, n}, the blocks Bk are 

Bk ^ a+ti{I + R) + (si + ak)I ~ PkR - IkU ■ 
Proof. From the definition of Bk and the computation of Ann, we have 

Bk = (51 + y)I - A„o + E AnA^'^''''^'^^ - I)- 

And, from the computation of v4„0j we obtain that Bk ~ a-^.fi{I + R) + sil + Dk, 
with 

n-l 

Dfc = ^A„,(e^-(™K-/). 
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Now, our problem has been reduced to calculate Dk- Using the explicit 
computation of Anj , we see that Dk satisfies 



n-l 



(2sin(jC/2))"+^ 
The coefficient of the sum can be written as 

a_(7 - e^i^'^'+'K) - a+R (e-^^'^ - e'^*=^) . 
Notice that, using the equalities 

e-UJ<) + e'^iJC) = 2/ cos jC, and 
^i(iki+j)c + g-j(ifc/+J)C = 2 [/cos jfcCcos iC + «^sin jfcCsin jQ] , 

we may cancel terms from the sum for j and n — j . In this way, we obtain 
that the matrix is 

a_ (/[I — cos fcjCcos jC] — zJ[sinjfcCsin jX]) — a+-R[cos jC ~ cos jfcC] 
h (2sin(iC/2))"+^ ■ 

Hence, we may write D/. as = a^I — — 7fei J with 

El — cos kj( cos j'C 

(2sin(,C/2))"+^' 
cos — cos j/cC 

(2sinO-C/2))"+^' 

Esin jfcC smj( 
(2sin(iC/2))"+- 

Finally, we conclude that ak, Pk and 7fe coincide with the definitions in the 
proposition from the equalities 

^ 1 - cos((fc + 1)jC) 

(2smOC/2)) ^ 
^ 1 - cos((fc - 1)jC) 
(2sm(jC/2)) ^ 
o ^ ^, sin^(jfcC/2)-sin^(jC/2) , - 

+ ^ (2sin(jC/2))"+' +^ ^ 



/3fe = Q!+ ^ 



Proposition 19 For fc e {l,n — 1}, we ftave that Bn-i = Bi and 

I n {si + fi + na-) - {^Y^"^ fia -(t)^^^A** 
Bi = — ^^y^"^ ^Q, si + ai + {a + l)fj, aii 
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Proof. Prom the proof of the previous proposition, we have that 

n 

A„jeJ(*^+-^)f =a+{I + R)n + (si + ai)/ - /3ii? - 7ii J. 
And, since 0i =0 and ai =71, then 



3{ii+J)i ^ ( Si + "1 + 2a+/x aii 
—aii si + ai 



Moreover, since Rvi = vi, then 

n^/^Ao„ui = -n^/^/i(a_f;i + a+vi) = ji 



-a 
i 



From the definition of Bi we get the result. Finally, using the computation of 
Bi, we may prove that RBiR — Bi, and then that i?„_i = RBiR = Bi. ■ 

Clearly, the sums Sk are positive and satisfy Sk = Sn+k = Sn-k- To analyze 
the bifurcation points, we need the following recursive formula for Sk- 

Lemma 20 Let Sk be defined as the sum Sk but with a — 2 instead of a. Then, 
the sums Sk satisfy the recurrence formulae 



Sk+i - Sfc = (2fc + l)si -"^Sh 



h=l 

Proof. We write the sum Sk as 

2«s - V ^ 1 ~ cos(A:jC) 

' ^4^sin«-i(jC/2) l-cos(iC) ■ 

Using geometric series, we have 

1 - cos(fciC) _ 1 - e'J'^g 1 - e-'i'"^ _ ^ 
1 - cosfjC) ~ l-e*J« l-e-*-Jf ~ 

Now, we may cancel common terms from Sk+i and Sk as 

n — 1 ^ A; 



" " 1 



Finally, since 

k 



^ e^i'^<^ J2 cosi/iC = (2fc + l)-4^sin20X/2), 

h= — k h=—k h=l 
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then 



k n—1 . 2/ 



Sk+i - Sk - {2k + - F F ^^^y^^ = {2k + 1)51 - y 



h=l 



The idea of using geometric series is taken from [T], where it is used to 
calculate Sk for the vortex case a = 1. Iterating this result we obtain the 
equalities 

Sk+i - 2sk + Sfe-i = 2si - Sk (10) 

and 

k-i / ; \ fc-i 

Sk = J2^si+i -si) = J2i {21 + l)si - ^ Sh = k^si - Isk-i. (11) 

1=0 1=0 \ h=l / 1=1 

4.1 General potential 

Now, we only need to find the bifurcation points, that is, the points where <7k{p-) 
changes sign for the general potential ([IJ. For k = n, we have that /3„ — — si 
and an = ot-Si, then ef i?„ei — {a + l)(/i + si) and 

o-„(^) = sgn{p. + si) 

Proposition 21 The sign a\ is 

ai{n) = sgn{hiii{ji + si){ji - pLi)), 

where fj,i = —ai/bi, with 

ai = (si + 2ai) (2si + na — n) and &i = (a + 1) (2si + 2ai — n) . 

Proof. We get the result from the fact that det Bi can be factored as follows 

2 det i?i , N ^ N / „ s /„ N 

— ^ fi{a + 1) (2si + 2ai - n) + {si + 2ai) (2si + na - n) 

^{fl+Sl) 

= 6i(/i-/xi). 



Remark 22 Notice that ai and an change sign at — Si, then 77i(— si) = and 
rjn{—si) — ±2. Nevertheless, there are two explicit bifurcations at —si one with 
symmetry Di and another one with D„. Indeed, the bifurcation with symmetry 
Di is made of the translations of a, (0 + r, e*'» + r, e'"^ + r), with oj — 0. The 
bifurcation with symmetry Dn is made of the homotheties of a, (0, e'^r, e™^r), 
with uj — 0. 

In addition, since a\ changes sign at fi — 0, then r]i{0) — ±2 . Therefore, 
there must be a bifurcation, with symmetry Di, at ii = 0. As the central body 
has mass zero, fi — 0, then the bifurcation has no physical meaning since it is 
made of the solutions {r, e*^, e™^), with fi — 0. 
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Proposition 23 For k E {2, [71/2]}, the signs (Jk are 



(^k{p) = sgn {fj, - ilk) 



where fik 



ak/bk, with 



7fc - /3fc and bk (a + 1) (si + + Pk)- 



Proof. For k = n/2, we have that 7^/2 = and 




,/2ei = si + a„/2 ~ P'nii + (a + 



For fc ^ {l,n/2,n}, the determmant of is 



From the definitions of cxk and we have that 



si + Ofe + ^fc = Sfc + (a_/2) (sfc+i + 2sfc + Sk-\ - 2si) . 



Using the equality (|10p and the fact that 4sfc— sj. is positive, we get the inequaUty 
Sfe+i + 2sfc + Sfc_i > 2si. Consequently, the factor 6^ is positive and we may 
conclude the result. ■ 

From the discussion in the previous remark, the true bifurcations are found 
at /ife for k e {1, [n/2]}. From the bifurcation theorem we have the following: 

Theorem 24 For k G {1, [n/2]}, let h be the maximum common divisor of 
k and n. If jik is different from —Si, and fij for the other j £ [l,n/2] H hN, 
then, from fj,k, there is a global bifurcation of relative equilibria with maximal 
symmetry Dh- 

By maximal symmetry Dh we mean that the local branch has symmetry Dh 
but not for a bigger group Dp. 

By global bifurcation we mean that, whenever the branch is admissible, the 
branch returns to other bifurcation points and the sum of the local degrees at 
these bifurcation points is zero. The branch is inadmissible when the parameter 
or the norm goes to infinity, or when the branch ends in a collision solution. 

Remark 25 Notice that these results are applicable only for n > 3, since the 
irreducible representations of the definition ^ are not consistent for n = 2. 
Nevertheless, the case n = 2 was analyzed in the same spirit in a previous 
remark. For instance, we have proved that there is a bifurcation of relative equi- 
libria with symmetry Di from fii = — (2a + Si) /(a + 1). Also we did calculate, 
for the vortex problem, that /ii = —5/4 and, for the body problem, fii = —17/12. 
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4.2 (n + l)-vortex potential 

Here we give a short description of the bifurcation points for the vortex problem, 
since, in this case, we can calculate explicitly the bifurcation points /ij.. 

Proposition 26 For a = 1, we have that 

Sk — k{n — k)/2. 

Proof. For a = 1, we have si — {n— l)/2. In addition, we may calculate Sk as 

n—1 n—1 

Sk - 2^ sin2(fcjC/2) = ^(1 - cos(fcjC)) = n. 
j=i i=i 

Therefore, from the formula (ITlT) . we have that 

fc-i 

Sk ^ {n ~ I) /2 - n'^l ^ k{n - k) /2. 

1=1 

■ 

From the definitions with a = 1, we have a_ = 0, a+ = 1, ak = 0, 7fe = 
and Pk = Sk — si- Since /ife — Sk/2 — si, for k e {2, [n/2]}, then 

^j, = (-fc2 + nfc - 2n + 2) /4. 

And, for fc = 1, we have 

Ail = = (n - 1)2/4. 

Consequently, the bifurcation point = ~l/2 is always negative and = 
(n — 7)/4 is positive only for n > 8. Moreover, since jik is increasing in n for 
fc > 3, then /ifc is always positive for fc > 4. Notice also that the bifurcation 
points /ifc are increasing in fc for fc G {2, [ri/2]}, and, as a consequence, the 
/ifc are different. 

Theorem 27 For n > 3, anc? eac/i fc G {1, [n/2]}, the polygonal relative 
equilibrium has a global bifurcation of relative equilibria from fXk with maximal 
symmetry Dh- 

The existence of the local bifurcation was proved before in the article [10] , 
with a normal form method. 

4.3 (ri + l)-body potential 

Notice that, for the (n+l)-body problem, the equations have a physical meaning 
only for /i > 0. Given that we cannot calculate explicitly the sums Sk in this case, 
we shall give an asymptotic computation of the sums and of the bifurcation 
points Ilk- 
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Proposition 28 For n big enough, the bifurcation point /ii is negative and fik 
is positive for k > 2. 

Proof. For the (n + l)-body problem a = 2. From the definitions, we have in 
this case a_ = 1/2, a+ = 3/2, ak = {sk+i + Sfe-i)/4, Pk = 3(sfc - Si)/2 and 
7fc = (■Sfc+i - Sfc-i)/4. 

Using integral estimates, it can be easily seen that si/n — >■ oo and that Sk/n 
is finite when n goes to infinity. Therefore, from the formula (|Tl1) . we have the 
limits Sfe/si — > A;^, when n goes to infinity. 

We have, for k > 2, that /Sk/si 3{k'^ - l)/2, ak/si (fc^ + l)/2 and 
7fc/si — >■ fc, when n — > oo. Therefore, from the definitions of Ok and bk, we 
obtain the limits bk/ Si — > 6fc^ and 

afc/s? = (1 + afc/si)2 - (fik/sif - (7fc/si)' ^ -fc' (2fc' - 5) . 

Consequently, the result follows from the fact that /ifc/si converges to the pos- 
itive hmit (2fc2 - 5) for k>2. 

For fc = 1, we have that ai/si — > 1, then we obtain the result from 

(si + 2ai)(2si+n) 
3 (2si + /q;i — n) 

■ 

In [1^, the bifurcation of the local branch from jik is proven for the {n + 1)- 
body problem. 

Remark 29 Given the numerical evidence of fXk, for instance see ]! ffl, it seems 
that /xi > for n G {3,4,5,6}, /i2 > for n > 10 and fik > for every 
A; > 3. The numerical evidence also suggests that the fj,k are increasing for 
k G {2, [ri/2]}. This is true at least in the limit when n 00, because 
{fJ.k+1 — Mfc)/si converges to the positive limit {2k + 1) /3. 

Theorem 30 Assuming the numerical evidence of the previous remark, from jii 
for n e {3, 4, 5, 6}, from 112 for n > 10, and from pik for each k £ {3, [?t./2]}, 
the polygonal relative equilibrium has a global bifurcation of relative equilibria 
with maximal symmetry Dh ■ 



5 dNLS 

The dNLS equations are 

iqj = h{\\qj\\^)qj + {q.j + i - 2q.j + q^^i), 

where qj € C represents the oscillator and h is the nonlinear potential. We wish 
to study a finite circular lattice, that is, a lattice of oscillators for j e {1, n}, 
with periodic conditions qj = qj+n- 

The solutions of the form qj = e"**Xj, with Xj constant, are called relative 
equilibria. In order to obtain the amplitude as a parameter, we need to change 
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coordinates, with qj = fie'^^^Xj. In this manner, we have that the values Xj form 
a relative equilibrium when 

—ujXj — h(\^Xj'^)xj + {xjj^i — 2xj + Xj^i). 

Remark 31 Given that the lattice is integrable for n — 1 and n — 2, we shall 
look for bifurcation of relative equilibria for n > 3. Actually, according to J^, it 
is possible to find all the bifurcation diagrams of the relative equilibria for n < 4. 
Notice that the relative equilibria are known as breathers when they are localized. 

The starting point is a relative equilibrium which looks like a rotating wave 
and is the equivalent of the polygonal relative equilibrium in the n-body prob- 
lem. We give next a condition which needs to be satisfied by the potential for 
the existence of this rotating wave. 

Proposition 32 Define aj ~ e*-''», with C, — 2Ti/n, then a = (ai,...,a„) is a 

relative equilibrium if 

uj = 4sin2(C/2) - h{fi^). 
Proof. Since aj+i — 2aj + aj^i — — 4sin^(C/2)aj, then 

V,^{a) = {lu + h{fj,^) - 4sin2(C/2))a,. 

■ 

Remark 33 Note that the existence of the rotating wave is determined by a 
non-homogeneous relation between the amplitude /i and the frequency w. This is 
different from the n-body problem, where the existence of the relative equilibrium 
is determined by a homogeneous relation. 

In order to show the similarities with the n-body problem we change to real 
coordinates. Let x = {xi,...,Xn) € K^" be the vector of positions, then the 
relative equilibria are critical points of the potential 

1 " 

where Xj = Xj+n and H{x^ /i) is a function such that WH{x) — ujx -f- h{\ij,x\^)x. 

From the point of view of the symmetries, there is practically no difference 
with the definitions of the n-body problem ([2]). The unique difference is in the 
fact that we are not including the coordinate xq of the n-body problem. So, in 
this case, the group Z)„ acts on R^" as p{j){xi, ...,Xn) = (a;-y(i), a;7(„)) and 
the group 0{2) = U kS^ in a similar way. The fact that the gradient is 
I?„-equivariant follows from the periodicity conditions Xj = Xj+n- Moreover, it 
is well known that the potential is invariant when we rotate the phases of all 
oscillators, so the gradient WV is 0(2)-equivariant. 
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As a consequence, we may adapt the results of sections two and three. Ac- 
tuaUy, as in Proposition ([7]), in this case the blocks are given by 

n 

Sfe = ^A„,e^(^'^-^+'^)^ (12) 

for k e {1, n/2, n}, and the signs a{fi) are defined as before in ([5]). Fur- 
thermore, since, in this case, there is no collision points, then the bifurcation 
is inadmissible only when the parameter /i or the norm of the branch goes to 
infinity. 

We wish to describe briefly the meaning of the symmetries Dh ^ for the 
dNLS equations. Due to Xj — e^^^'^'^^'^^ Xj_,_n/h- then the solutions look like 
rotating waves composed of h identical waves, each one formed by n = n/h 
oscillators which satisfy the reflection symmetry xj = Xn-j- An example of 
relative equilibria with symmetry Dh is 

Xj = (l + £sin2j(7r/n))e'(2'^/''^ 

Given that most of the work is already done, we shall focus our attention on 
finding the bifurcation points. 



5.1 General potential 

Again, the first step is to find the submatrices Aij of D^V{x) at aj. 

Proposition 34 The submatrices Anj are Anj — I for j E {1, n — 1}, Anj — 
for j ^ {l,n — l,n} and 

Ann = (-2cosC)/ + 2/i2/i'(^^)dmg(l,0). 

Proof. As the coupling is linear and only between adjacent oscillators, then 
Anj ^lioij e {l,n-l}, Anj ==Ofor j ^ and A„„ ^ D'^H{an)-2I. 

Let xq — {x,y), since \JH{xq) = ujxq + /i(|/ia;o|^)xo, then 

D^H{xo) = {u + h{\^,x,\^))I + 2^?h\\^,x,\') ^J;^ ^! 

Since a is an equilibrium when uj + h{fi'^) — 4sin^(C/2), then at a„ = (1,0) we 
have 

D^H{an) ^ 4sin^{C/2)I + 2n^h'{\n\^)diag{l,0). 

Hence, we conclude the result from the equality 4sin^(C/2) — 2 = — 2cosC. ■ 
Now we may calculate the blocks Bk from (fT2|) . 

Proposition 35 Define and -fk as 

ak = 4cosCsin^ k(/2 and 7^. ~ 2sinfcCsinC. 
Then, the blocks Bk are 

Bk = -aki + jk{iJ) + 2^i^h'{\^l\^)diag{l, 0). 
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Proof. Using the explicit computation of Anj , we have 

Bk = (^e'^'''^+-^^^ + 6-^'''^+''^'^^ + i~2cosO I + 2n^h'{fi^)diagil,0). 
Then, from the equahties 2cosC(cosfcC — 1) = — afe and 

g(»fc/+j)c ^g-(«fe/+./)C ^ (2cosA:CcosC)/+ (2sinfcCsinOjJ, 

we obtain the form oi B^- ■ 

Now, it remains to find the bifurcation points. Since a„ = 0, then we have 
eiBnBi — 2fi'^h'{p'^) and 

cr„ — sgn h'{p^). 

Since, for n = 3, we have ai = —71 < 0, then cti = sgn h'{fi^). As for n = A, 
we have ak = and 71 7^ 0, then cri — —1 and a2 — sgn Given that, 

in our examples, h'{y?) does not change sign, there are no bifurcation points 
for n = 3, 4. 

Consequently, we shall focus our attention only on the cases n > 5, for 
k S {1, n/2}, where we can assume 

ttfe = 4cosCsin^ fcC/2 > 0. 

Proposition 36 For n > 5 and k e {l,...,n/2}, the sign ak can change sign 
only for the solutions of ^'^h' {p'^) — 5k, with 

Sk = {al - 7fe)/(2afe). 
Moreover, we have Si < 0, ^2 = and Sk > for k G [3, n/2] D N. 
Proof. For k G [l,n/2) n N, we have 

detBk^al~Jl-2ak^l^h'{^i^) 
= 2ak {Sk - li^h'ii?)) . 

Since ak > 0, then ak = sgn [Sk — (x^h'i^^^y). 

When k = n/2, we have 7„/2 = 0, then i?„/2ei = — a„/2 + 2^'^h'{p?). 
Thus, <T„/2 = sff"- {iJ-^h'lii^) — an/2) changes sign only for the solutions of 
fi'^h'ifj.'^) = a„/2/2 = (5„/2. 

Finally, since Sk has the sign of 

al - 7fc = 16(sin2 fcC/2 - sin^ C) sin^ K/^, 

then Sk has the sign of sin^ k(/2 — sin^ (. ■ 

From the bifurcation theorems (fT4|) and (|15p . we have the following result:. 

Theorem 37 For each simple solution of ^J^k^'{^J^k) ~ ^k, from the amplitude 
jJLk we have a global bifurcation of relative equilibria with symmetry Dh, where 
h is the maximum common divisor of k and n. 

Remark 38 Actually, we may analyze more complex lattices whenever we pre- 
serve the symmetries. For instance, we may consider nonlinear coupling and 
coupling with distant oscillators. 

Now, we wish to give two typical examples. 
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5.2 The Schrodinger cubic potential 

For the cubic Schrodinger potential, we need to set h{x) = x. In this case 
h'{ii^) = 1 and a„ = 1. 

Then, for n > 5, the sign crfc(/L<) changes only when /Xfe = \/5k, if 5k is positive. 

As wc have proven before that 5k is positive when n > 6, for k G {3, [n/2]}, 
and, since the numbers 5k are increasing in k, then the fjtk are increasing for 
e {3, [n/2]}. 

Theorem 39 For the cubic Schrodinger potential, for n G {6,7,..}, for each 
k e {3,..., [n/2]} there is a global bifurcation of relative equilibria with maximal 
symmetry Dh from the amplitude \fSk. 

5.3 A saturable potential 

For a saturable potential, we need to set /i = (1 + x)~^. In this case, h^ji^) = 
-(1 + M^)"^, cTn = -1 and 

mW) = -m^(i + /^V 

is a function with range (—1/4,0) and a single minimum at /ti^ = 1. Therefore, 
there are two zeros, /z_ G (0, 1) and G (1, oo), of the equation = 5k, 

when 5k G (-1/4,0). 

Since we have proven before that ^fe > for fc > 2, there is no bifurcation 
for fc > 2 and it remains only to analyze the case k = 1. Since 

6i = 2{sm^ C/2 - sin^ ()/ cos C -5> 0" 

when n — >■ oo, then 5i G (—1/4,0) if n is big enough. Indeed, we obtain 
numerically that G (—1/4,0) for n > 16. 

Theorem 40 For the lattice with saturable potential, for n > 16, from the am- 
plitudes /x_ G (0,1) and 11+ G (l,oo) there is a bifurcation of relative equilibria 
with maximal symmetry Di . 
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